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CHAPTER  I 


INTRODUCTION 


1.1  Description  of  Draft  Tube  Surge 


During  the  past  decade,  considerable  attention  has  been  directed 
toward  the  production  of  energy  by  pollution-free  methods.  As  a result 
of  this  interest,  emphasis  has  been  placed  on  the  development  of  high 
head  pump-turbine  units  because  of  their  production  of  relatively 
inexpensive  energy.  Pump-turbines  in  general  are  used  to  satisfy  peak 
power  requirements  using  the  pump-storage  concept.  Because  of  the  high 
variability  of  required  power  production  of  a peaking  power  plant,  pump- 
turbines  are  required  to  run  over  a wide  range  of  operating  conditions. 

A phenomenon  called  draft  tube  surge  is  associated  with  hydroelectric 
turbines  and  pump-turbines  which  operate  at  off  design  conditions.  Draft 
tube  surge  is  a hydrodynamic  instability  which  occurs  as  a result  of 
excess  rotation  remaining  in  the  water  as  it  leaves  the  turbine  runner 
and  enters  the  draft  tube.  Figure  1 is  a schematic  of  a turbine  cross 
section.  Due  to  the  excess  rotation,  a swirling  flow  (axial  motion  and 
rotational  motion  at  the  same  time)  is  established  in  the  draft  tube  of 
the  turbine.  When  the  flux  of  angular  momentum  entering  the  draft  tube 
becomes  sufficiently  large  compared  to  the  flux  of  linear  momentum  in  the 
draft  tube,  a reversal  of  flow  direction  occurs  along  the  draft  tube 
centerline.  In  the  field  of  fluid  machinery  this  condition  is  known  as 
inception  of  draft  tube  surge.  A rope-like  helical  vortex  develops 
around  the  region  of  reversed  flow  and  precesses  about  the  draft  tube 
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Schematic  of  a turbine  cross  section 
illustrating  the  geometry  of  the  flow 
passage  in  the  region  of  the  wicket 


centerline.  Cassidy  and  Falvey  ]_  2_/  reported  that  the  presence  of  the 
helical  vortex  produced  unsteady  pressures  in  the  draft  tube.  Falvey  and 
Cassidy  j_  4 _/  later  observed  the  frequency  of  the  unsteady  pressure  to  be 
very  regular,  indicating  that  the  helical  vortex  precessed  about  the  tube 
centerline  at  a reasonably  regular  rate.  Along  with  this  periodic  vibra- 
tion in  the  draft  tube,  the  surges  are  responsible  for  many  other 
undesirable  operating  characteristics  in  hydroelectric  power  plants. 

These  undesirable  characteristics,  which  may  occur  individually  or  in 
combination,  are  noise,  vibrations,  variation  in  power  output,  vertical 
movement  of  the  turbine  runner  and  shaft,  and  pressure  pulsations  in  the 
penstock. 

The  draft  tube  surge  phenomenon  is  complicated  by  many  factors, 
including  vaporization  of  water,  geometry  of  the  turbine  and  draft  tube, 
and  the  dynamic  characteristics  of  the  penstock  and  electrical  network. 

The  vortex  core  which  arises  as  a result  of  the  swirling  flow  in  the 
draft  tube  has  been  observed  to  be  filled  with  air,  water,  or  water  vapor. 
Falvey  J_  5_/  noted  that  the  fluid  phase  of  the  vortex  core  was  related  to 
the  tailwater  elevation  of  the  turbine.  For  high  tailwater  elevations, 
the  vortex  core  was  observed  to  be  filled  with  water,  whereas  for  the 
lower  tailwater  elevations,  water  vapor  and  sometimes  air  was  found 
present  in  the  vortex  core.  The  two-phase  flow  which  occurred  for  low 
tailwater  conditions  was  observed  to  produce  pressure  surges  which  were 
somewhat  reduced  in  magnitude. 

The  geometric  shape  of  the  draft  tube  has  also  been  observed  to  be 
a governing  influence  on  surging  characteristics.  Palde  t_  9 _/  conducted 
an  extensive  study  on  the  influence  of  draft  tube  shape  on  surging 
characteristics  and  found  that,  for  a given  turbine  installation,  the 
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geometry  of  the  draft  tube  could  affect  the  range  of  surging,  the 
frequency  band  of  the  surge,  and  the  amplitude  of  the  resulting  pressure 
pulsations. 

In  an  effort  to  coordinate  the  studies  of  draft  tube  surge, 
attempts  were  made  to  develop  predictive  terms  of  surge.  Murakami  J_  8 _/ 
noted  that  the  frequency  and  force  of  surge  vibrations  could  be  predicted 
when  the  dimensions  of  the  draft  tube  and  turbine  runner  and  load  condi- 
tions of  the  turbine  were  known.  Cassidy  and  Falvey  j_  2_/  later  showed 
that  frequency  and  pressure  pulsations  of  surge  could  be  correlated  to 
the  amount  of  angular  momentum  in  the  draft  tube.  In  a later  paper  by 
Falvey  and  Cassidy  j_  4_/,  dimensionless  parameters  were  developed  to 
describe  surge.  It  was  assumed  that,  for  a particular  draft  tube  shape, 
the  root-mean-square  amplitude  |apJ  and  frequency  (f)  of  the  unsteady 
pressure  surge  are  both  functions  of  the  fluid  density  (p)  and  kinematic 
viscosity  (v),  draft  tube  diameter  (D)  and  length  (L) , discharge  (Q)  and 
flux  of  angular  momentum  (ft) . Standard  techniques  of  dimensional 
analysis  gave  the  following  functional  relationships: 


D4AP  ftp  L 

_2  * *1  2 ’ D * 

pQ  PQ 


and 


(1) 


(2) 


2 

where  Rg  is  the  Reynolds  number  4Q/ttDv.  Inspection  of  ftD/pQ  shows  that 
this  represents  the  ratio  of  angular  momentum  flux  to  linear  momentum 
flux  in  the  draft  tube  and  is  defined  as  the  momentum  parameter.  The 
ratio  D AP/pQ  is  called  the  pressure  parameter,  and  fD  /Q  is  called  the 
frequency  parameter.  Experimental  studies  proved  that  a high  correlation 


i 


*r  «*  « 4 


exists  between  all  three  parameters.  For  any  hydraulic  turbine  installa- 
tion, relations  of  pressure  parameter  versus  momentum  parameter, 
frequency  parameter  versus  momentum  parameter,  and  frequency  parameter 
versus  pressure  parameter  can  be  developed,  provided  that  a method  is 
available  to  determine  the  momentum  parameter.  Palde  j_  10 _/  developed  a 
graphical  technique  called  the  swirl  momentum  method  to  predict  the 
momentum  parameter  and  the  occurrence  of  draft  tube  surge  for  a given 
wicket  gate  geometry.  A more  accurate  means  of  predicting  momentum 
parameter  with  respect  to  wicket  gate  opening  was  developed  by  Yocum 
f_  17 / through  the  use  of  a computer  model. 

1.2  Prior  Attempts  to  Eliminate  Draft  Tube  Surge 

Because  of  the  detrimental  effect  draft  tube  surge  has  on  hydrau- 
lic power  plants,  various  attempts  have  been  made  over  the  years  to 
alleviate  the  problem,  some  dating  back  to  1880  and  earlier. 

Probably  the  most  frequently  mentioned  method  used  to  reduce  draft 
tube  surge  is  the  admission  of  air  into  the  runner  or  draft  tube  (Figure 
2).  Both  location  and  manner  of  air  injection  have  been  found  to  be 
important.  Air  has  been  admitted  upstream  of  the  runner,  between  the 
runner  and  the  wicket  gates,  through  the  runner  core,  and  through  the 
draft  tube  wall.  For  low  tailwater  elevation  installations,  the  pressure 
near  the  draft  tube  centerline  is  less  than  atmospheric  and  all  that  is 
necessary  for  air  admission  is  proper  ventilation  to  the  vortex  core. 
However,  for  higher  tailwater  levels,  the  use  of  an  air  compressor  is 
necessary  to  inject  air  into  the  draft  tube.  Consequently,  air  admission 
may  be  costly  because  the  most  severe  surging  occurs  at  high  tailwater 
elevations.  Air  admission  was  found  to  reduce  the  pressure  fluctuation 
in  the  draft  tube  of  several  installations.  However,  other  installations 
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using  air  admission  have  reported  an  increase  in  surge  j_  5_7.  This  is 
due  to  the  dependence  of  surge  on  cavitation  number  which  was  observed  by 
Gerich  and  Raabe  [_  6_/. 

Because  surge  is  caused  by  excess  rotation  remaining  in  the  flow 
after  leaving  the  runner,  attempts  have  been  made  to  eliminate  or  reduce 
the  excess  rotation  by  placing  fins  parallel  to  the  draft  tube  axis  and 
attached  to  either  the  draft  tube  wall  or  the  runner  core  (Figure  3A). 

If  the  fins  extend  far  enough  across  the  draft  tube  throat  to  touch  each 
other,  they  are  know  as  flow  splitters  (Figure  3B)  /_  5_/.  Fin  geometries 
varying  the  number  of  fins  along  with  length  and  depth  of  fins  have  been 
tried,  but  all  with  limited  success.  The  disadvantages  of  fins  and  flow 
splitters  are  loss  in  efficiency,  necessity  to  guard  the  appentures 
against  cavitation  damage,  high  frequency  noise,  and  structural  problems 
involved  in  installation.  Because  of  these  disadvantages,  the  use  of 
fins  is  normally  restricted  to  small  units  ]_  5_/. 

Other  investigations  with  model  test  facilities  have  shown  that  a 
hollow  cylinder  which  resembles  a barrel  with  the  ends  cut  out  and  placed 
concentrically  with  the  draft  tube  will  reduce  the  amplitude  of  the  pres- 
sure surges  (Figure  3C) . The  coaxial  hollow  cylinder  has  been  noted  to 
reduce  pressure  pulsations  up  to  70%.  By  placing  the  hollow  cylinder  in 
the  draft  tube,  the  vortex  core  is  coaxially  contained  in  the  draft  tube. 
Falvey  j_  5_/  noted  that  a vortex  core  which  forms  coaxially  with  the 
draft  tube  produces  no  surge.  However,  a major  problem  with  the  coaxial 
hollow  cylinder  is  structural.  Supports  which  are  necessary  to  keep  the 
cylinder  centered  in  the  draft  tube  would  have  to  be  strong  enough  to 
resist  large  lateral  forces,  yet  small  enough  to  minimize  flow  blockage. 
The  problem  of  cavitation  of  the  appendages  would  also  have  to  be 


considered. 


FLOW  FINS  FLOW  SPUTTERS  COAXIAL  HOLLOW  CYLINDER 

(A)  (B)  (C) 

Figure  3.  Appendages  placed  in  the  draft  tube  in  an 
effort  to  contain  or  reduce  the  rotational 


Probably  the  oldest  approach  used  in  an  attempt  to  reduce  draft 
tube  surge  is  to  replace  the  vortex  core  present  in  the  draft  tube  with  a 
solid  body  of  rotation.  The  solid  body  can  be  either  an  integral  part  of 
the  draft  tube  as  in  the  case  of  the  Moody  Spreading  Cone  (Figure  4A)  or 
attached  to  the  runner  core  (Figure  4B)  j_  5_/.  The  reasoning  behind  the 
use  of  extensions  to  the  runner  core  is  that  the  elimination  of  the 
reverse  flow  area  reduces  unsteadiness  by  streamlining  the  overall  flow. 
However,  the  use  of  runner  core  extensions  allows  unfavorable  velocity 
head  recovery  in  the  draft  tube  and  produces  high  lateral  forces  on  the 
runner  extension.  In  addition  to  these  disadvantages,  an  adjustable 
diameter  extension  would  be  necessary  for  optimal  results  since  the 
diameter  of  the  reverse  flow  region  varies  from  0.2  to  0.8  draft  tube 
diameters. 

Of  all  the  above  methods  to  reduce  draft  tube  surge,  only  limited 
success  has  been  obtained  with  any  of  them.  Each  of  the  methods,  to  some 
degree,  causes  efficiency  loss,  results  in  cavitation  damage,  induces 
structural  problems,  or  requires  expensive  auxiliary  equipment. 


effort  to  streamline  draft  tube  flow. 


CHAPTER  II 


THE  REDUCTION  OF  DRAFT  TUBE  SURGE  BY  INJECTION  OF 
COUNTER  ROTATIONAL  FLUID  IN  THE  DRAFT  TUBE 
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2.1  Description  of  the  Fluid  Injection  Method 

The  subject  of  this  thesis  is  directly  related  to  a newly  proposed 
method  to  eliminate  draft  tube  surge  in  hydroelectric  turbines  and 
pump- turbines.  The  method  considered  here  consists  of  a series  of  flush- 
mounted  nozzles  in  the  draft  tube  placed  immediately  downstream  of  the 
turbine  runner.  The  nozzles  will  inject  high  energy  fluid  in  a periph- 
eral direction  with  respect  to  the  draft  tube  axis  and  into  the  draft 
tube  counter  to  the  rotating  motion  of  the  swirling  flow.  The  intent 
behind  this  method  is  momentum  exchange  between  the  injected  fluid  and 
the  swirling  fluid  exiting  from  the  turbine  runner.  With  the  proper 
amount  of  fluid  injection,  swirl  in  the  draft  tube  could  be  reduced  below 
the  point  of  draft  tube  surge  inception.  A schematic  of  the  proposed 
arrangement  is  shown  in  Figure  5 where  the  high  energy  injection  fluid  is 
bypassed  from  the  spiral  casing  at  the  turbine  inlet.  Although  not 
indicated  in  the  figure,  two  separate  rows  of  nozzles  could  be  provided 
with  one  row  of  nozzles  capable  of  injecting  fluid  in  a direction 
opposite  to  that  of  the  other.  The  reduction  of  positive  or  negative 
swirl,  which  is  encountered  at  below  design  and  above  design  operating 
conditions  respectively,  would  be  possible  with  this  type  of  arrangement. 
An  appropriate  system  of  valving  activated  on  the  basis  of  wicket  gate 
opening  could  control  the  flow  to  the  nozzles  and  control  the  number  of 
nozzles  discharging. 


' 
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2.2  Desirable  Characteristics  of  Injection  of  Counter  Rotational  Fluid 
The  reduction  of  draft  tube  surge  by  fluid  injection  has  several 
advantages  over  the  other  methods  which  have  been  used  to  reduce  surge. 

In  the  proposed  method,  the  nozzles  will  be  flush-mounted  on  the  wall  of 
the  draft  tube  and,  therefore,  a minimum  amount  of  flow  disturbance  will 
be  present  when  the  use  of  the  nozzles  is  not  required.  The  presence  of 
the  nozzles  should  not  cause  any  loss  in  efficiency  or  cavitation  resist- 
ance of  the  machine  when  there  is  no  injection,  provided  that  cavitation 
does  not  occur  in  the  slots  when  the  jets  are  not  in  use.  As  the  turbine 
operates  progressively  further  away  from  design  conditions,  the  swirl 
nozzles  could  be  activated  to  provide  an  increasing  amount  of  counter- 
swirl to  the  draft  tube.  By  using  such  a system,  an  infinite  number  of 
off  design  conditions  which  are  encountered  during  turbine  shutdown  and 
start-up  could  be  handled.  A second  row  of  nozzles  capable  of  injecting 
fluid  opposite  to  that  of  the  other  row  would  permit  the  reduction  of 
negative  swirl.  Also,  this  system  will  not  require  the  placement  of 
appendages  in  the  draft  tube  that  will  be  subject  to  fluctuating  hydro- 
dynamic  loads.  The  feasibility  of  this  proposal,  however,  hinges  upon 
the  amount  of  fluid  which  must  be  bypassed  from  the  spiral  casing  to  the 
injection  nozzles  in  order  to  eliminate  a given  surging  condition.  Any 
fluid  which  is  bypassed  represents  an  energy  loss  across  the  machine  and 
therefore  produces  an  efficiency  loss.  Therefore,  an  evaluation  of  the 
amount  of  fluid  which  must  be  drawn  from  the  high  pressure  side  of  the 
turbine  and  discharged  into  the  draft  tube  is  of  primary  importance.  An 
estimate  of  the  ratio  of  bypassed  flow  to  flow  through  the  turbine  can  be 
obtained  on  the  basis  of  momentum  considerations. 
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2.3  Prediction  of  Bypass  Flow  Requirements 

For  any  turbine  installation,  the  amount  of  angular  momentum 

present  in  the  draft  tube  can  be  expressed  by  the  dimensionless  momentum 
2 

parameter  ftD/pQ  for  any  given  operating  condition.  The  basic  expression 
for  power  can  be  shown  as 

P - wT  , (3) 

where  u is  the  angular  velocity  of  the  turbine  runner,  and  T is  the 
torque  of  the  turbine  runner.  The  expression 

T - - 02  , (4) 

where  £2^  - £2^  is  the  rate  change  of  moment  of  momentum  of  the  fluid  as  it 
passes  through  the  turbine  runner  can  be  substituted  into  Equation  (3)  to 
obtain  the  expression 

1 * * S • (5) 

This  equation  can  be  nondimensionalized  by  multiplying  both  sides 

2 

by  D/pQ  . Rearranging  terms.  Equation  (5)  takes  the  form  of  the  equation 
obtained  by  Palde  j_  13 _/ : 

V _ V . _PD_ 

2 2 2 * 

pQ  pQ  pmQ 

This  expression  shows  very  simply  that  the  amount  of  angular  momentum 
present  in  the  draft  tube  for  any  given  condition  is  equal  to  the  amount 
of  angular  momentum  induced  by  the  turbine  wicket  gates  minus  the  torque 


imparted  to  the  turbine  runner  and  converted  to  mechanical  energy.  To 
maximize  the  effect  of  fluid  injection  on  surge  reduction,  the  amount  of 
bypass  fluid  required  to  eliminate  surge  must  be  minimized.  Therefore, 


the  momentum  parameter  of  a given  surging  condition  must  be  reduced  to 
the  critical  value  of  momentum  parameter  which  describes  inception  of 
draft  tube  surge. 
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The  momentum  parameter  or  swirl  parameter  shall  be  defined  as  S, 

where 


S 


(7) 


The  relative  influence  of  the  various  operating  characteristics  of  the 
turbine  with  respect  to  swirl  parameter  becomes  more  apparent  if  the 
swirl  parameter  can  be  expressed  as  a ratio  of  angular  velocity  to  axial 
velocity  in  the  draft  tube.  Figure  6 illustrates  the  nomenclature  and 
stations  considered.  From  Figure  6,  the  draft  tube  flow  (Q)  and  flux  of 
angular  momentum  (J2)  can  be  determined  at  the  exit  of  the  wicket  gates 
(WG)  as 


/ v' 

J WG 


Q =>  / V dA  = V 2irr.B 

1 R R l 


and 


0 - pf  VR  ( riV9 . ) <**  ■ pQrlY 

WG  1 1 


(8) 


(9) 


where  B is  the  wicket  gate  depth, 
velocity  at  the  wicket  gate  exit. 
Equation  (7)  yields: 

pQriYD  riVD 
S - 2s - 

PQ  Q 


and  V is  the  rotational  component  of 
0i 

Substitution  of  Equation  (9)  into 


(10) 


where  D is  the  draft  tube  diameter.  The  presence  of  a given  amount  of 
swirl  in  the  draft  tube  above  the  point  of  inception  of  surge  will 
generate  a pressure  pulsation  of  some  prescribed  frequency.  The  swirl 
in  the  draft  tube  shall  be  generated  by  a set  of  wicket  gates  as  indica- 
ted by  Fiyire  6.  Conservation  of  momentum  is  assumed  between  the  wicket 
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gates  and  the  draft  tube  throat.  From  Figure  6,  the  quantity  V 

nr 

represents  the  mass  averaged  axial  velocity  in  the  draft  tube.  The 

assumption  is  made  that  the  radial  distribution  of  tangential  velocity  at 

a given  section  in  the  draft  tube  is  represented  by  a potential  vortex. 

On  this  basis,  the  mass  averaged  quantity  V , at  a given  section  of  the 

°T 

draft  tube,  would  occur  at  the  diameter  D associated  with  the  mean  area 
of  the  draft  tube.  This  can  be  expressed  in  algebraic  form  as 


r v « — v = — — v 
1 9i  2 9t  2yr  9i 


(id 


If  the  results  of  Equation  (11)  are  substituted  into  Equation  (10),  the 
relation  between  swirl  and  velocity  ratio  is  found  to  be 


D2  V 


e„ 


2~\J~2  V 


ttD 


or 


V7 


* V 


1 

*T 


(12) 


To  evaluate  the  amount  of  bypass  fluid  required  to  eliminate  a 
particular  surge  condition,  it  shall  be  assumed  that  the  swirl  parameter 
present  in  the  draft  tube  before  injection  will  be  reduced  to  some  criti- 
cal swirl  parameter  ( S ] (*  demotes  critical  condition) . The  quantity 

of  bleed  flow  required  to  obtain  the  critical  swirl  parameter  ( S ] in 
the  draft  tube  can  be  evaluated  by  equating  the  difference  between  the 
angular  momentum  from  the  turbine  discharge  and  that  of  the  bypass  fluid 
to  the  angular  momentum  at  the  critical  condition.  On  this  basis,  the 


relation 
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(13) 


QTV9  ~ QBVB  “ QT  V9 
T T 

is  obtained,  where  the  subscript  T refers  to  the  turbine  throat,  and 
subscript  B refers  to  the  bypass  fluid  flow.  Equation  (13)  can  be  non- 

"ft 

dimensionalized  by  Q„  V to  give: 
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Ve. 


0 v 

, VB  B 


0 V 
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Q V 


S* 


VT 


(14) 


The  velocity  of  the  bleed  fluid  V can  be  expressed  as  V2gH>  where  H is 

B 

the  turbine  head.  Using  the  continuity  principle  and  Equation  (12), 

Equation  (14)  can  be  rearranged  to  provide  a ratio  of  bypass  flow  to 

* 

draft  tube  flow  in  terms  of  the  swirl  parameters  S and  S . The  predic- 
tive equation  becomes: 


2.22  VA  * 
*T 


B 


A 
— * 


- S 


— * 


(15) 


T 2gH 

Examining  the  above  equation  for  a given  swirl  parameter  (S),  it 
is  seen  that  the  numerator  of  the  above  expression  will  be  essentially 
constant,  independent  of  turbine  head.  This  is  based  on  the  assumption, 
for  cavitation  purposes,  that  the  velocity  in  the  draft  tube  throat  at 
the  design  condition  for  conventional  turbine  installations  has  some 
upper  limit.  Therefore,  Equation  (15)  indicates  that  the  required  ratio 
Q„/Q_  to  eliminate  a particular  swirl  parameter  (S)  decreases  as  some 
reciprocal  function  as  head  increases.  Since  the  ratio  Q„/Q_  represents 

O 1 

a loss  in  efficiency,  this  would  indicate  that  the  decrease  in  efficiency, 
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CHAPTER  III 

INJECTION  NOZZLE  GEOMETRY  STUDIES 

3.1  Optimization  of  Nozzle  Geometry  for  Effective  Reduction  of  Draft 
Tube  Surge 

The  successful  reduction  of  draft  tube  surge  by  fluid  injection 
depends  upon  the  efficient  exchange  of  momentum  between  the  injection 
nozzle  flow  and  the  rotational  flow  of  the  swirling  condition.  It  would 
appear  that  the  geometry  of  the  injection  nozzle  would  have  an  effect  on 
the  efficient  exchange  of  momentum.  Therefore  an  extensive  series  of 
tests  were  performed  to  determine  the  most  effective  nozzle  geometry  for 
fluid  injection.  Various  geometries  with  different  length- to-width 
ratios  along  with  circular  geometries  were  tested.  The  effect  of  fluid 
injection  angle  with  respect  to  a tangent  to  the  draft  tube  wall  was  also 
investigated.  The  results  of  these  studies  indicate  that  the  optimum 
nozzle  geometry  is  dependent  on  the  head  and  flow  characteristics  of  a 
given  pump-turbine  installation.  Table  I (page  37)  is  a summary  of  the 
nozzle  geometries  tested. 


3. 2 Test  Facility  Description 

The  investigations  for  this  thesis  were  conducted  in  an  air  flow 
facility  as  shown  in  Figure  7.  Studies  conducted  by  Cassidy  ]_  1_/  and 
Palde  J_  9_/  using  a similar  facility  have  shown  that  the  use  of  air  as  a 
working  fluid  has  been  quite  successful  in  characterizing  the  swirl 
parameter  in  draft  tube  flow  and  predicting  the  occurrence  of  draft  tube 
surge.  The  main  air  supply  was  provided  by  a variable  speed  centrifugal 
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Figure  7.  Schematic  of  test  facility  used  in  the  draft  tube 
surge  study. 
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blower  and  was  measured  by  an  orifice  meter  arrangement.  The  metered  air 
was  diffused  as  it  entered  a stilling  chamber  where  screens  were  used  to 

reduce  flow  turbulence  across  the  stilling  chamber  flow  field.  The  air 

entered  the  cylindrical  draft  tube  and  nozzle  test  section  radially 

through  wicket  gate  type  swirl  vanes.  The  angle  between  the  vanes  and  a 

radial  line  could  be  set  at  any  value  between  0°  (radial  position)  and 

82.5°  (closed  position).  Later  in  this  thesis,  it  will  be  shown  that  the 

air  flow  angle  through  the  wicket  gates  is  directly  proportional  to  the 

momentum  parameter.  The  nozzle  test  section  provided  three  points  of 

fluid  injection  equally  spaced  in  a circumferential  plane  perpendicular 

to  the  draft  tube  centerline.  Figure  8 is  a photograph  of  the  nozzle 

test  section  with  nozzle  //5  in  place.  Injection  fluid  provided  by  an 

auxiliary  centrifugal  blower  and  measured  by  an  orifice  meter  arrangement 

entered  the  manifold  of  the  injection  nozzle  test  section  and  passed 

through  the  nozzle  openings  into  the  draft  tube.  Any  ratio  of  injection 

fluid  flow  (bypass  flow)  (Q  ) to  fluid  flow  through  the  wicket  gates 

5 

(turbine  flow)  (Q^)  could  be  provided  by  the  proper  setting  of  a butterfly 
valve  connected  to  the  auxiliary  air  supply  blower.  After  passing  through 
the  nozzle  test  section,  the  air  flow  continued  through  the  cylindrical 
draft  tube  (L/D  = 4.56).  Pressure  taps  and  probe  holes  were  located  at 
various  locations  along  the  length  of  the  tube.  The  studies  were  con- 
ducted with  a cylindrical  draft  tube  and  also  with  an  elbow  type 
(Fontenelle)  draft  tube. 

3.3  Test  Facility  Instrumentation 

The  unsteady  pressure  produced  by  the  swirling  flow  in  the  draft 


tube  was  monitored  at  the  last  pressure  tap  (L/D  ■ 4.40)  on  the  draft  tube 
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by  a dynamically  calibrated  1 psi  (0.689N/cm  ) differential  pressure 
transducer.  The  dynamic  calibration  of  the  pressure  transducer  was 
performed  by  Dr.  Frank  Archibald  of  the  Applied  Research  Laboratory  of 
The  Pennsylvania  State  University  as  described  by  Dr.  Archibald  in 
Appendix  A.  The  differential  pressure  transducer  was  connected  to  the 
draft  tube  pressure  tap  with  a 4"  (10  cm)  long  piece  of  surgical  tubing 
having  a 1/4"  (6.35  mm)  inside  diameter  and  was  referenced  to  atmospheric 
pressure.  To  obtain  the  magnitude  of  the  pressure  surge,  the  transducer 
signal  was  sent  to  a DISA  (Type  55  D35)  RMS  meter  through  a Krohn-Hite 
(Model  3342)  line  filter.  The  line  filter  had  a characteristic  48  dB  per 
octave  rolloff  and  was  set  to  pass  everything  between  20  and  120  Hz  for 
all  tests.  The  frequency  of  the  pressure  surge  was  obtained  by  passing 
the  raw  signal  to  a Spectral  Dynamics  Real  Time  Analyzer  (Model  SD-301) 
which  was  operated  in  conjunction  with  a Spectral  Dynamics  Ensemble 
Averaging  Unit  (Model  SD-302).  Knowing  the  pressure  surge  RMS  value  and 
the  frequency  of  the  surge,  both  pressure  parameter  and  frequency  param- 
eter could  be  calculated. 

Velocity  profiles  in  the  draft  tube  and  behind  the  wicket  gates 
were  obtained  by  the  use  of  two  three-hole  pressure  probes.  A 1/8"  (3.18 
mm)  diameter  prism  probe  was  used  through  a 3/16"  (4.76  mm)  wide  arc  slot 
in  the  wicket  gate  assembly  for  traversing  across  the  trailing  edge  of 
the  wicket  gates  to  determine  the  velocity  profile  and  fluid  flow  angle. 
Figure  9 shows  typical  flow  angle  variation  behind  the  wicket  gates  with 
respect  to  circumferential  location.  Similarly,  Figure  10  shows  typical 
variations  of  velocity  with  respect  to  circumferential  location.  In  the 
draft  tube,  a 1/4"  (6.35  mm)  diameter  wedge  probe  was  used  to  investigate 
velocity  profiles  and  pressure  distributions  across  the  diameter  of  the 


CIRCUMFERENTIAL  LOCATION  (deg) 


Figure  9.  Variation  of  flow  angle  behind  the  wicket 
gates  with  respect  to  circumferential 
location. 
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draft  tube.  The  probe  was  inserted  into  the  draft  tube  through  probe 

holes  provided  at  different  locations  along  the  draft  tube  centerline. 

Both  probes  were  used  in  conjunction  with  a multi-channel  scanner  valve 

2 

and  a 0.5  psi  (0.345  N/cm  ) differential  pressure  transducer.  Figure  11 
is  a schematic  of  the  draft  tube  instrumentation  used  in  the  study. 


3.4  Preliminary  Data  Required  Before  Testing 

Before  any  significant  data  could  be  taken,  the  wicket  gate 
assembly  had  to  be  calibrated  to  produce  a relation  between  momentum 
parameter  and  wicket  gate  angle.  Instead  of  using  wicket  gate  angle,  an 
arbitrary  linear  scale  was  used  to  determine  gate  opening.  When  the 
wicket  gates  were  closed,  the  scale  read  0.0.  When  the  gates  were 
radial,  the  scale  read  9.5.  Therefore,  a relation  between  momentum 
parameter  and  wicket  gate  setting  could  be  developed. 

To  simplify  the  calibration  process,  the  relation  for  swirl  param- 
eter as  described  by  Equation  (10)  was  examined.  The  angle  between  a 
radial  line  and  the  velocity  vector  V shown  in  Figure  6 was  defined  to  be 
0,  and  therefore,  from  Figure  6: 

VR  = V cos  9 (16) 

and 

V.  =■  V sin  0 . (17) 

9i 

Using  Equation  (8),  the  relation  for  volume  flow  rate,  along  with 
Equations  (16)  and  (17)  and  substituting  into  Equation  (10),  it  could  be 
shown  that 
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Figure  11.  Schematic  of  draft  tube  instrumentation. 
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For  the  test  facility,  the  values  of  D and  B indicated  by  Figure  6 were 
known  constants.  Equation  (18),  therefore,  could  be  reduced  to 

S = Kx  tan  9 , (19) 

where  ■ 0.557.  Equation  (19)  shows  that  momentum  parameter  (swirl 
parameter)  is  only  a function  of  fluid  flow  angle  at  the  wicket  gate 
trailing  edge  for  the  test  facility.  To  calibrate  the  wicket  gate 
assembly,  a series  of  measurements  were  taken  across  the  wicket  gates 
trailing  edge  at  various  wicket  gate  settings.  The  arithmetic  mean  flow 
angle  was  determined  across  five  wicket  gates  (75  data  points)  for  each 
wicket  gate  setting.  With  this  information,  a relation  between  swirl 
parameter  and  wicket  gate  setting  was  developed  and  is  shown  in  Figure  12. 
It  was  also  necessary  to  determine  the  relations  of  pressure  parameter  and 
frequency  parameter  as  functions  of  momentum  parameter  for  the  particular 
test  facility.  Palde  ]_  9_/  reported  the  effect  of  draft  tube  shape  upon 
surging  characteristics  and,  therefore,  a separate  set  of  relations  had  to 
be  developed  for  both  the  cylindrical  and  elbow  type  draft  tubes.  For 
various  momentum  parameters,  the  RMS  value  and  frequency  of  the  pressure 
surge  were  determined  by  the  instrumentation  described  earlier.  From  this 
data,  the  necessary  relations  could  be  developed.  Since  the  test  facility 
used  for  this  study  was  very  similar  to  that  used  by  Palde  j_  9_/,  it  was 
desirable  to  get  similar  relations  of  pressure  parameter  and  frequency 
parameter  as  a function  of  momentum  parameter.  Figures  13  and  14  show  the 
relations  obtained  in  this  study  for  the  cylindrical  and  elbow  type  draft 
tubes,  respectively.  On  the  same  figures,  Palde's  data  are  plotted  to 
show  that  the  data  are  in  fairly  good  agreement  with  the  previously  deter- 
mined relations. 
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Figure  14.  Comparison  of  the  surging  characteristics 
of  the  elbow  type  draft  tube  with  Palde's 
data.  Circles  show  present  test  results. 
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Preliminary  investigations  were  also  performed  in  the  draft  tube 
flow  field.  Velocity  profiles  were  developed  for  a series  of  momentum 
parameters.  Both  the  axial  velocity  component  and  the  rotational 
velocity  component,  nondimensionalized  by  the  mass  average  axial  velocity 
in  the  draft  tube,  were  plotted  as  functions  of  nondimens ional  radius. 
Figure  15  shows  the  characteristic  curves  for  swirl  parameters  of  0.24, 
0.41,  0.80,  and  1.18.  It  can  be  seen  in  the  axial  velocity  distribution 
that  a region  of  reverse  flow  occurs  for  all  swirl  parameters  and  becomes 
increasingly  larger  as  swirl  parameter  increases  in  magnitude.  The  axial 
velocity  distributions  also  approach  a potential  vortex  at  the  outer 
region  and  the  inner  region  takes  the  form  of  a solid  body.  The  rota- 
tional velocity  distribution  shows  that  the  product  V r increases  for 
increasing  swirl  parameter.  At  S = 0.24,  there  is  no  surge;  yet,  there 
is  considerable  rotation  in  the  flow.  The  inception  of  surge  is  very 
close  to  S * 0.33. 

3.5  Nozzle  Geometry  Optimization  Procedure 

The  optimization  of  nozzle  geometry  was  performed  in  two  parts. 
First,  the  effects  of  angle  of  fluid  injection  (<j>)  with  respect  to  a 
tangent  line  to  the  draft  tube  wall  in  a circumferential  plane  was  inves- 
tigated. The  nozzle  shape  and  area  were  held  constant  while  various 
injection  angles  were  tested  for  effectiveness  in  reducing  draft  tube 
surge  for  S = 0.80  and  S = 1.18.  The  values  of  S * 0.80  and  S ■ 1.18  were 
selected  for  testing  because  these  are  values  of  swirl  parameter  which  are 
commonly  encountered  in  field  installations.  From  these  data,  an  optimum 
angle  of  injection  was  selected  and  the  remaining  series  of  tests  used 
this  optimum  angle  of  injection  which  investigated  the  effects  of  various 
nozzle  geometry  shapes.  The  effects  of  length-to-width  ratio  along  with 
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varying  nozzle  area  were  investigated  for  their  effectiveness  in  the 
reduction  of  draft  tube  surge.  All  nozzle  optimization  tests  for  injec- 
tion angle  effects  and  nozzle  shape  effects  were  performed  with  the 
cylindrical  draft  tube  attached  to  the  test  facility. 
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CHAPTER  IV 

EXPERIMENTAL  RESULTS 

4.1  Optimization  of  Angle  of  Injection 

The  first  step  in  the  nozzle  geometry  optimization  procedure  was 
to  determine  the  most  effective  angle  of  fluid  injection  for  the  reduc- 
tion of  draft  tube  surge.  For  this  series  of  tests,  the  nozzle  area  and 
shape  were  held  constant  while  the  angle  of  injection  (<(>)  was  studied  at 
four  different  values,  namely,  0°,  15°,  30°,  and  45°.  Table  I defines 
the  angle  of  injection  with  respect  to  the  draft  tube  and  also  lists  the 
various  nozzle  geometries  investigated  in  this  study.  Nozzle  #5  was  used 
for  the  optimization  of  the  fluid  injection  angle.  It  should  also  be 
pointed  out  that  nozzle  area  A^^  in  Table  I is  the  sum  of  all  three  nozzle 
areas  injecting  into  the  draft  tube.  Also,  A^,  is  defined  as  the  area  of 
the  draft  tube  throat. 

To  determine  the  optimum  angle  of  injection,  each  angle  was  studied 
to  see  how  it  affected  the  characteristics  of  the  surging  condition.  Of 
utmost  importance  is  the  effect  of  angle  of  injection  on  the  reduction  of 
pressure  parameter.  The  amount  of  bypass  fluid  required  to  eliminate 
surge  at  swirl  parameters  of  0.80  and  1.18  was  measured  for  each  angle  of 
injection.  The  Spectral  Dynamics  Real  Time  Analyser  was  used  to  determine 
when  surge  was  gone.  Figure  16  shows  a typical  series  of  plots  received 
from  the  analyser.  Since  dBs  are  related  to  volts  by  the  expression 

Volts  = -20  Lob  dB  (20) 


and  volts  are  related  directly  to  pressure  through  the  transducer 
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constant,  it  can  be  seen  in  Figure  16  that  surge  pressure  was  systemati- 
cally reduced  by  increasing  amounts  of  injection  fluid.  Notice  that  at 
an  injection  ratio  Q /Q  of  16%,  surge  disappeared  from  the  analyser 
spectrum.  Figure  17  shows  the  effect  of  increasing  fluid  injection  upon 
the  magnitude  of  the  unsteady  pressure  in  the  draft  tube  for  all  four 
injection  angles  at  both  S * 1.18  and  0.80.  The  plots  show  that  up  to  an 
angle  of  45°,  there  is  no  significant  difference  in  reduction  of  pressure 
parameter  for  any  given  angle  of  injection.  At  S = 1.18,  the  flow  ratio 
required  for  the  elimination  of  surge  is  essentially  constant  at  16%  for 
all  injection  angles.  A similar  statement  can  be  made  for  S » 0.80. 

Velocity  profiles  were  also  obtained  in  the  draft  tube  at  probe 
hole  L/D  = 4.40  for  each  angle  of  injection  at  the  point  of  surge  elimina- 
tion. Figure  18  shows  that  for  the  four  angles  of  injection,  the 
rotational  velocity  profile  and  the  axial  velocity  profile  are  essentially 
the  same  at  these  conditions.  With  this  evidence,  it  was  concluded  that 
for  an  angle  of  injection  up  to  45°,  there  was  no  measurable  difference 
in  the  effective  reduction  of  draft  tube  surge  using  nozzle  geometry  //5. 
For  this  study,  it  was  assumed  that  this  was  true  for  all  nozzle 
geometries.  However,  in  order  to  continue  the  experiments  on  nozzle 
geometry  shape,  an  angle  of  injection  had  to  be  selected.  An  injection 
angle  of  30°  was  used  for  the  remaining  series  of  tests. 

4.2  Optimization  of  Nozzle  Shape 

Table  I is  a list  of  the  nozzle  geometries  tested.  The  first  nine 
nozzle  shapes  in  Table  I were  initially  studied  for  area  ratio  effects. 
However,  nozzles  it 2,  it 5,  and  itl  comprised  a series  of  nozzles  where  the 
area  ratio  was  constant  and  the  shape  varied.  Nozzle  shapes  10  through  15 
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were  used  to  compliment  the  study  of  length-to-width  effects  upon 
constant  area  ratio  nozzles. 

Each  nozzle  was  placed  into  the  test  facility  and  tested  for  its 
effectiveness  in  the  reduction  of  surge  for  swirl  parameters  of  0.80  and 
1.18.  A nozzle  which  required  less  bypass  fluid  to  eliminate  surge  was 
declared  more  efficient  than  the  other.  A progressive  look  at  the  reduc- 
tion of  surge  in  terms  of  pressure  parameter  with  respect  to  increasing 
flow  ratio  was  studied  for  each  nozzle.  For  all  cases,  an  increase  in 
the  amount  of  injected  fluid  would  result  in  a decrease  in  the  magnitude 
of  the  pressure  parameter  and  frequency  parameter  as  earlier  illustrated 
by  Figure  16.  Figure  19  shows  the  effect  of  increasing  flow  ratio  on 
pressure  parameter  for  swirl  conditions  of  1.18  and  0.80.  The  nozzle 
shapes  used  in  these  plots  were  selected  to  show  the  effect  of  decreasing 
area  ratio  on  the  amount  of  bypass  fluid  required  to  eliminate  a given 
surging  condition.  The  last  data  point  on  each  curve  defines  the  exact 
amount  of  fluid  required  to  eliminate  surge  for  the  given  nozzle  shape. 


Figure  19  shows  that  for  a given  ratio  of  Q_/Q_.,  the  reduction  in  the 

D 1 

surge  pressure  amplitude  is  greater  as  A^/A^,  decreases.  This  results 
directly  from  the  fact  that  to  inject  a given  ratio  of  QD/Q_  into  the 
draft  tube,  the  turbine  head  and,  hence,  the  nozzle  fluid  velocity  must 


increase  as  A^/A^  decreases.  Since  the  momentum  of  the  injected  fluid 
increases  as  the  square  of  the  spouting  velocity  from  the  nozzles,  smaller 
flow  ratios  are  required  to  eliminate  surge  as  the  head  of  the  turbine 
increases.  This  was  suggested  by  Equation  (15). 

The  effect  of  various  nozzle  length-to-width  ratios  on  surge 


reduction  was  studied  for  area  ratios  of  1.27,  2.54,  and  5.09,  all 
expressed  in  percent.  For  a swirl  parameter  of  1.18  and  a constant  area 


ratio  of  2.54,  the  amount  of  bypass  fluid  required  for  surge  elimination 
varied  from  16.1%  to  19.8%  with  respect  to  the  various  length- to-width 
ratios  tested.  However,  at  A^/A,j,  = 1.27,  the  amount  of  bypass  fluid  only 
varied  from  12.5%  to  14.1%,  and  at  A^/AT  = 5.09,  the  variation  was  even 
smaller  (21.7%  to  22.6%).  Figure  20  is  a plot  of  flow  ratio  versus  area 
ratio,  which  includes  the  results  of  all  the  nozzle  shapes  tested  for  S = 
1.18.  Figure  21  is  a similar  figure  only  with  S = 0.80.  Figures  20  and 
21  show  that  for  a constant  area  ratio  of  2.54%,  there  was  an  observable 
difference  in  the  amount  of  injection  fluid  required  to  eliminate  surge 
for  various  L/W  ratios  at  S = 1.18  and  to  a lesser  degree  at  S * 0.80, 
However,  at  constant  area  ratios  of  1.27%  and  5.09%  for  both  S = 1.18  and 
0.80,  the  amount  of  required  injection  fluid  measured  for  the  various  L/W 
ratios  was  considered  to  be  within  the  limits  of  experimental  error.  It 
was  further  noticed  that  only  nozzle  //7  in  Figure  20  indicated  that  vari- 
ability in  length-to-width  of  a nozzle  would  cause  a significant  increase 
or  decrease  in  effective  reduction  of  surge.  Since  that  data  point  was 
only  one  in  twenty  examined  for  L/W  effects,  it  was  concluded  that  no 
appreciable  difference  was  observed  in  the  ratio  of  Q_ / Q_  required  to 
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eliminate  surge  at  a given  A^/A^,  for  the  various  nozzle  geometries  tested. 
This  result  suggests  that  the  choice  of  nozzle  geometry  to  be  used  for  a 
given  Ajj/A^  to  eliminate  draft  tube  surge  would  consist  of  selecting  that 
geometry  which  would  produce  minimum  hydraulic  losses  and  permit  the 
easiest  fabrication.  On  this  basis,  a circular  nozzle  would  seem 
appropriate. 

4.3  Effects  of  Injection  on  Surge  Characteristics 

Figures  16  and  17  were  presented  to  show  the  effects  of  injection 
on  pressure  parameter  and  frequency  parameter  of  a given  surging 


A^/At  (percent) 


Figure  20.  Injection  flow  required  to  eliminate  surge 

for  all  nozzles  tested  using  the  cylindrical 
draft  tube,  S = 1.18. 
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Figure  23.  Mixing  of  injected  flow  with  draft  tube  flow  with 
respect  to  draft  tube  length. 
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CHAPTER  V 

COMPARISON  OF  TEST  RESULTS  WITH  PREDICTIONS 

5.1  Development  of  Empirical  Relations 

A relation  to  predict  the  ratio  of  bypass  flow  to  turbine  flow  was 
previously  derived  and  presented  in  Equation  (15) . This  relation  was 
developed  assuming  ideal  momentum  transfer  between  the  injected  fluid  and 
the  rotational  flow  in  the  draft  tube.  An  empirical  relation  applicable 
to  pump- turbines  can  also  be  derived  using  the  experimental  data  pre- 
sented in  Figures  16  and  17.  The  bypass  flow  can  be  expressed  as  the 
product  of  the  mass  averaged  velocity  through  the  nozzle  times  the  nozzle 
area  or 


qB  - K2(28H)1/2Ali 


(21) 


where  is  a nozzle  coefficient  that  is  assumed  to  have  the  value  of 


0.9.  The  flow  rate  through  the  turbine  can  be  defined  by  the  familiar 
power  expression  as 
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(22) 


By  combining  Equations  (21)  and  (22),  an  expression  for  flow  ratio 
develops  into 

l3/2 
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(23) 


where  y * 62.4  lb/ft^  (1000  Kg/M^),  and  g = 32.2  ft/sec^  (9.81  m/sec^). 
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Multiplying  both  numerator  and  denominator  by  A^,  and  substituting  in  the 
known  values  of  given  quantities.  Equation  (23)  becomes: 


or 


Q TJ  ^fr  0/0 

» 450.7  (English) 

V [3987  "^^“3/2]  <“■“«> 


(24) 


With  the  experimental  data  shown  in  Figures  16  and  17,  empirical 
relations  to  predict  Qg/Q^  as  a function  of  A^/A^,  were  developed  for  S * 
1.18  and  0.80  and  are  plotted  on  the  figures.  The  solid  data  points  were 
selected  for  the  simple  linear  regression  of  the  two  variables  because 
these  points  represented  nozzles  that  were  circular  or  had  an  L/W  ratio 
equal  or  close  to  1.0.  The  other  data  points  were  not  used  in  order  to 
eliminate  the  effects  of  L/W  ratio  which  were  not  detected  in  this  study 
but  could  possibly  still  exist.  For  S * 1.18,  a least  squares  fit  of  the 
data  produced  the  relation 
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and  similarly,  for  S = 0.80: 

0.366 


B 


0.473  -r- 


(26) 


Equations  (25)  and  (26)  provide  a convenient  means  of  eliminating 
from  the  right-hand  side  of  Equation  (24)  for  their  respective 
swirl  conditions.  Therefore,  if  a swirl  condition  of  S = 1.18  is  con- 
sidered, Equation  (25)  can  be  substituted  into  Equation  (24)  to  give: 
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In  a similar  fashion,  the  substitution  of  Equation  (26)  into  Equation 
(24)  provides  a relation  of  Qtt/Q  for  a swirl  parameter  of  0.80, 
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Observation  of  Equations  (27)  and  (28)  shows  that,  given  the 
geometric  and  performance  data  of  a turbine  installation,  the  required 
amount  of  bypass  fluid  to  eliminate  surge  can  be  predicted  for  the  surge 
conditions  of  S = 1.18  and  0.80. 

Data  similar  to  that  shown  in  Figures  20  and  21  were  also  obtained 
with  an  elbow  type  draft  tube  in  place  of  the  cylindrical  draft  tube. 

The  elbow  type  draft  tube  tested  was  similar  to  the  Fontenelle  configura- 
tion described  by  Palde  /_  9_/.  Empirical  relations  were  also  developed 
for  this  data  set  and  are  shown  in  Figure  24,  plotted  with  the  data  for 
S = 1.18  and  0.80.  Using  these  empirical  equations  for  the  elbow  type 
draft  tube  and  Equation  (24),  relations  similar  to  Equations  (27)  and 
(28)  can  be  developed  in  likewise  manner  for  the  elbow  type  draft  tube. 
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5.2  Estimates  of  Bypass  Flow  for  a Turbine  Installation 

In  the  report  by  Palde  J_  10 _/,  geometric  and  performance  data  were 
given  for  the  model  of  the  Grand  Coulee  Third  Power  Plant  Units.  There- 
fore, the  Grand  Coulee  Third  Power  Plant  Units  will  be  used  as  an  example 
for  the  calculation  of  the  bypass  fluid  required  to  eliminate  surge. 

With  the  geometric  and  performance  data,  which  are  shown  in  Figure  25, 
the  required  bypass  flow  rates  to  eliminate  surge  in  the  turbine  model 
for  S ■ 1.18  and  0.80  can  now  be  determined  using  Figure  21  and  Equations 
(27)  and  (28),  respectively.  If  we  assume  a speed  coefficient  $2  of 
0.80,  then  the  necessary  performance  data  can  be  determined  from  Figure 
25.  ""able  II  is  a summary  of  the  calculations  necessary  to  predict 

Q /Q  for  swirl  conditions  of  S = 1.18  and  0.80  using  the  empirical  data. 

B T 

Calculations  shown  are  for  the  cylindrical  draft  tube,  although  the 
results  for  the  elbow  draft  tube  are  shown  also.  Table  III  shows  the 
calculations  required  for  solution  of  the  theoretical  relation.  Equation 
(15) , and  Table  IV  is  a comparative  summary  of  the  results  of  Tables  II 
and  III. 

Table  IV  Indicates  that  the  empirically  derived  quantity  of  flow 
ratio  required  to  eliminate  surge  at  a swirl  parameter  of  S * 1.18  is 
about  2.5  times  greater  than  that  predicted  by  assuming  ideal  momentum 
transfer. 

The  differences  between  the  estimated  empirical  quantity  and  the 
estimated  ideal  quantity  of  Q /Q  for  S **  0.80  is  smaller.  The  discrep- 
ancy  between  the  ideal  results  and  the  empirical  results  is  most  likely 
due  to  poor  assumptions  made  in  the  development  of  Equation  (15).  The 
assumption  that  ideal  momentum  transfer  would  occur  in  the  test  facility 
draft  tube  between  the  injection  flow  and  rotational  flow  is  probably  not 
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1.18 

0.80 


0.1319 

0.1129 
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valid.  Therefore,  the  results  of  the  empricial  data  should  not  be 
expected  to  match  with  the  ideal  data.  Table  IV  also  shows  that  the 
shape  of  the  draft  tube  affects  the  efficiency  of  the  injection  process. 
More  importantly,  the  rotational  velocity  profile  in  the  draft  tube  was 
assumed  to  take  the  form  of  a potential  vortex  and  the  axial  velocity 
distribution  was  assumed  to  be  uniform.  Figure  15  shows  that  the  assumed 
velocity  distributions  were  poor  assumptions. 


CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

The  results  of  the  experimental  studies  of  injection  nozzle 
geometries  to  implement  the  reduction  of  draft  tube  surge  by  fluid 
injection  have  led  to  the  following  observations: 

1.  Injection  of  high  energy  fluid  counter  to  the  rotational  flow 
in  the  draft  tube  will  reduce  draft  tube  surge. 

2.  As  the  ratio  of  injected  fluid  Q /Q  is  increased,  both  the 
surge  frequency  and  the  magnitude  of  the  surge  amplitude  decreases. 

3.  An  angle  of  fluid  injection  up  to  45°  does  not  affect  the 
efficiency  of  the  fluid  injection  scheme. 

4.  The  use  of  various  nozzle  geometries  with  a fixed  area  ratio 
indicates  no  measurable  change  in  the  effectiveness  of  the  nozzle  with 
respect  to  eliminating  surge. 

5.  For  field  installation,  the  selection  of  nozzle  geometry  shape 
would  be  primarily  based  on  hydraulic  efficiency  and  ease  of  fabrication. 

6.  Empirical  data  are  presented  which  permit  estimating  the  area 
of  injection  nozzles  and  the  quantity  of  bleed  fluid  required  for  turbines 
of  specified  performance  characteristics. 

7.  For  a given  swirl  parameter,  the  bypass  flow  required  to 
eliminate  surge  decreases  as  higher  head  turbines  or  lower  specific  speed 


units  are  considered. 


— ■ 


CHAPTER  VII 

RECOMMENDATIONS 

With  the  completion  of  this  research  on  draft  tube  reduction, 
there  emerges  some  areas  of  interest  which  require  additional  attention. 
The  following  is  a list  of  recommendations  for  future  investigation: 

1.  Appendix  B shows  the  results  of  preliminary  tests  concerned 
with  the  effect  of  test  apparatus  geometry  on  the  pressure  amplitude  of 
draft  tube  surge.  The  short  study  shows  that  a stilling  chamber  placed 
upstream  or  downstream  of  the  wicket  gate  assembly  has  in  most  cases 
caused  a reduction  in  surge  pressure  amplitude.  Therefore,  stilling 
chambers  placed  upstream  or  downstream  of  the  turbine  runner  of  a field 
installation  may  be  an  effective  means  of  reducing  draft  tube  surge.  The 
geometry  of  the  stilling  chamber  and  its  location  relative  to  the  turbine 
runner  are  probably  important  considerations.  Further  study  is  required 
in  this  area. 

2.  The  theoretical  Equation  (15)  derived  in  section  2.3  to 
predict  the  amount  of  bypass  fluid  required  to  eliminate  a given  surging 
condition  has  been  proven  to  be  a poor  predicting  equation  by  the  results 
of  section  5.2.  Another  approach  to  developing  a predictive  equation  of 
required  bypass  rates  should  be  investigated  using  more  valid  assumptions. 

3.  Table  IV  shows  that  the  amount  of  bypass  fluid  required  to 
eliminate  a given  surging  condition  is  dependent  upon  the  shape  of  the 
draft  tube  (the  elbow  and  cylindrical  shape  are  compared  in  Table  IV). 
Further  investigations  are  required  to  determine  the  influence  of  draft 
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During  the  course  of  the  Draft  Tube  Surge  Investigation,  the 
question  was  raised  about  the  dynamic  sensitivity  of  the  transducer  being 
used.  The  static  sensitivity  of  this  transducer  was  found  to  be  0.0234 
psi/volt,  linear  over  the  dynamic  operating  range.  In  order  to  investi- 
gate whether  the  dynamic  sensitivity  was  different  from  the  static 
sensitivity  for  frequencies  up  to  100  Hz,  the  transducer  was  calibrated 
twice,  against  a condenser  microphone  in  both  cases. 

1 

The  frequency  response  of  each  Bruel  and  Kerr  condenser  microphone 
(1/2"  and  1/4")  was  flat  from  'vO  to  10  kHz,  as  evident  from  the  individual 
calibration  curves  provided  by  the  manufacturer.  The  sensitivity  of  each 
microphone  (including  amplifier  system)  was  determined  by  using  a piston 
phone  calibrator. 

In  the  two  calibration  situations  to  be  described,  a Spectral 
Dynamics  Real  Time  Analyser  (Model  SD-301)  was  used  in  conjunction  with  a 
Spectral  Dynamics  Ensemble  Averaging  Unit  (Model  SD-302) . The  procedure 
was  first  to  spectrally  analyse  the  250  Hz,  124  dB  calibration  signal  from 
the  microphone  as  produced  by  the  piston  phone  calibrator  to  give  an 
equivalent  electrical  spectral  level.  Subsequent  spectral  analyses  of  the 
microphone  signal  used  this  tone  level  from  the  calibration  spectrum  and 
accounted  for  analyser  attenuator  settings  and  system  gains.  Spectra  of 
the  subject  transducer  were  calibrated  using  the  inherent  calibration  of 
the  spectrum  analyser;  i.e.,  input  attenuator  setting  0 dB  reference  is 
100  mv.  As  this  same  analyser  is  used  for  all  the  experimental  spectral 
data  analysis,  the  machine  internal  calibration  will  be  the  same. 

However,  the  analyser  electrical  calibration  has  been  checked  by  applying 
a known  (as  measured  by  a high  quality  voltmeter)  RMS  level  sine  wave.  In 


the  calibration  situations  to  be  described,  the  transducer  was  exposed  to 
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Che  same  pressure  field  as  Che  microphone,  and  using  Che  specCral 
analyser,  Che  eleccrical  oucpuc  of  che  cransducer  was  deCermined  over  Che 
frequency  range  of  inceresC.  The  Cransducer  sensiCiviCy  can  Chen  be 
decermlned;  i.e.,  S - P/psi//Volc. 

The  Cwo  calibradon  slcuaclons  used  were  a caviCy  of  uniform 
(spaclally)  unsceady  pressure  and  an  unsCeady  pressure  from  a Capping  in 
Che  Drafc  Tube  Surge  ApparaCus.  Figure  26  is  a diagram  of  che  pressure 
cavicy  showing  Che  microphone  and  Che  Cransducer  locacions  on  che  bo c tom 
of  Che  cavicy  and  che  loudspeaker  on  Che  Cop.  Since  Che  highesc  fre- 
quency of  inCerest  is  100  Hz,  which  corresponds  to  an  approximate  11-ft 
wave  length,  the  cavity  dimensions  are  much  less  than  a wave  length  and 
the  pressure  will  be  uniform  within  the  cavity.  This  presumes  that  the 
loudspeaker  cone  acts  as  a piston,  an  assumption  that  again  ought  to  be 
satisfied  for  this  small  diameter  speaker  over  the  frequency  range  of 
interest.  A white  noise  electrical  signal  (0-20  kHz)  was  applied  to  the 
loudspeaker  at  two  levels,  5 and  0.5  Vrms.  This  resulted  in  an  almost 
exact  20  dB  difference  in  the  pressure  level  over  the  0-100  Hz  range  as 
measured  by  the  microphone.  It  is  noted  that  the  lowest  level  was 
approximately  20  dB  overall  (linear  weighting  2-40  kHz)  above  the  noise 
background  of  the  main  bay  area. 

The  spectral  level  of  the  pressure  transducer  signal  also  showed 
the  same  approximately  20  dB  change  over  the  frequency  range  between  the 
two  input  white  noise  power  levels,  thus  confirming  response  linearity 
with  signal  amplitude.  Figure  27  compares  spectra  taken  with  the  1/2" 
microphone  and  pressure  transducer.  The  frequency  response  of  the 
pressure  transducer  appears  quite  flat,  i.e.,  little  difference  between 
its  spectral  response  and  the  microphone,  known  to  be  flat.  Note  that 
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Figure  26.  Schematic  diagram  of  constant  pressure  calibration 
chamber  showing  location  of  microphone,  pressure 
tapping,  and  loudspeaker. 
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although  the  pressure  in  the  cavity  (as  measured  by  the  microphone)  is 
not  uniform  with  frequency,  the  deviations  are  approximately  duplicated 
by  the  signal  from  the  pressure  transducer.  It  is  possible  to  derive  the 
frequency  response  of  the  pressure  transducer  by  finding  the  difference 
in  its  response  relative  to  the  response  measured  by  the  microphone.  For 
example,  to  find  the  transducer  sensitivity,  take  the  values  at  30  Hz. 

The  mean  sound  pressure  level  is  52  dB  and  the  electrical  level  of  the 
pressure  transducer  at  this  frequency  is  0.0001  Vrms;  therefore,  S = 

0.0117  psi/volt. 

In  order  to  dispel  any  doubt  that  the  pressure  transduced  can 
accurately  measure  the  pressure  from  a tapping  in  the  research  apparatus, 
a 1/4"  microphone  was  connected  with  an  approximately  2"  long  piece  of 
1/4"  diameter  plastic  tubing.  Figure  28  compares  the  spectra  of  the 
microphone  and  pressure  transducer  for  the  same  operating  condition,  i.e., 
surge  was  strongly  present.  The  spectra  were  overlayed  to  match  the  level 
at  the  top  of  the  dominant  "hat  stack."  Using  the  values  at  this  fre- 
quency Cv-30  Hz),  the  sensitivity  of  the  pressure  transducer  was  computed 
to  be  0.0265  psi/volt.  This  is  in  good  agreement  with  the  static  sensi- 
tivity. In  addition,  it  is  seen  that  the  pressure  transducer  signal 
agrees  quite  well  with  the  signal  from  the  microphone  from  20  to  80  Hz. 
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The  test  facility  for  the  draft  tube  surge  study  as  shown  in 
Figure  7 was  not  the  original  test  apparatus  considered  in  the  program. 

The  original  draft  tube  surge  test  facility  was  designed,  fabricated,  and 
assembled  as  shown  in  Figure  29.  The  main  difference  between  the  original 
test  facility  and  the  one  finally  used  for  testing  is  the  location  of  the 
air  supply  flow  fan.  In  the  original  test  facility  (Figure  29),  air  was 
drawn  from  the  open  atmosphere  into  the  test  section  and  dumped  into  a 
stillin'  chamber  by  means  of  a flow  fan  downstream  of  the  draft  tube.  In 
the  final  test  facility  (Figure  7),  air  was  blown  into  a stilling  chamber 
first,  then  passed  through  the  test  section  and  dumped  to  atmosphere. 

Examination  of  the  sketches  of  the  two  facilities  would  seem  to 
indicate  that  the  important  flow  characteristics  (pressure  parameter  and 
frequency  parameter)  produced  in  either  test  facility  would  be  essen- 
tially the  same.  It  was  considered  important  for  the  test  facility  to 
produce  relations  between  the  surge  parameters  and  the  momentum  parameter 
similar  to  those  obtained  by  Palde  ]_  13 _/.  Preliminary  tests  in  the 
original  test  apparatus  provided  a relation  between  frequency  parameter 
and  momentum  parameter  much  like  that  obtained  by  Palde.  However, 
pressure  parameter  as  a function  of  momentum  parameter  was  essentially 
constant,  which  radically  disagreed  with  Palde1 s data.  In  an  attempt  to 
obtain  pressure  parameter  data  similar  to  Palde’ s,  various  pressure 
pick-up  devices  were  tried  at  various  locations  along  the  draft  tube. 

After  repeated  failure  to  obtain  the  desired  data  in  the  original  test 
facility,  a test  facility  similar  to  that  used  by  Palde  (i.e.,  pushing 
air  through  the  draft  tube  and  dumping  to  atmosphere)  was  constructed  as 
shown  in  Figure  7.  Surging  characteristics  obtained  in  the  modified  test 
facility  agreed  very  well  with  Palde’s  data.  The  altered  test  facility 
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solved  the  problem  but  raised  a question.  Why  was  the  pressure  parameter 
much  lower  for  the  original  test  facility  compared  to  that  obtained  in 
Palde's  or  the  final  test  facility?  Velocity  profiles  obtained  in  the 
draft  tube  at  L/D  * 4.40  for  each  facility  at  various  momentum  parameters 
were  compared.  Figure  30  shows  that  the  velocity  profiles  are  quite 
similar  for  both  test  facilities. 

In  an  effort  to  answer  this  question,  a brief  series  of  tests  were 
performed  to  investigate  the  effects  on  pressure  parameter  of  dumping  the 
draft  tube  flow  into  a stilling  chamber  and  then  dumping  to  atmosphere. 
Tests  were  performed  on  both  the  cylindrical  and  elbow  type  draft  tubes. 
Tables  B1  and  B2  show  the  different  geometries  that  were  investigated. 

With  any  given  facility  geometry  tested,  a survey  of  pressure  parameter 
as  a function  of  momentum  parameter  was  obtained.  The  flow  rate  through 
the  test  facility  was  measured  by  an  orifice  meter  and  the  surge  pressure 
was  measured  through  a dynamically  calibrated  pressure  transducer  by  an 
RMS  meter  as  earlier  described  in  this  report. 

Figure  31  shows  the  results  of  this  investigation  for  the  cylin- 
drical draft  tube  and  likewise  Figure  32  corresponds  to  the  Fontenelle 
(elbow  type)  draft  tube.  The  numbers  on  the  curves  in  Figures  31  and  32 
correspond  to  the  geometries  so  numbered  in  Tables  B1  and  B2.  Comparison 
of  curves  (1)  and  (2)  in  Figure  31  shows  the  effects  of  a stilling  chamber 
below  the  draft  tube  on  pressure  parameter.  Looking  at  all  the  curves  of 
both  Figures  31  and  32,  it  is  obvious  that  the  use  of  the  stilling  chamber 
in  all  cases  reduced  the  pressure  parameter  magnitude.  Notice  that 
screens  placed  across  the  stilling  chamber  / curve  (5) , Figure  3^/  forced 
the  pressure  parameter  relation  of  the  original  test  facility  to  become 
almost  identical  to  those  of  Palde  and  curve  (2). 
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TABLE  B1 

TEST  APPARATUS  GEOMETRIES  INVESTIGATED  WITH  RESPECT  TO  SURGING 
CHARACTERISTICS  USING  THE  CYLINDRICAL  DRAFT  TUBE 


1.  STILLING  CHAMBER  AND  AIR  SUPPLY  DOWNSTREAM  OF  DRAFT  TUBE 


OPEN 

ATMOSPHERE 


WICKET 

GATE 

ASSEMBLY 


DRAFT  TUBE 


AIR 

SUPPLY 


•STILLING  CHAMBER 


2.  NO  STILLING  CHAMBER,  UPSTREAM  AIR  SUPPLY 


AIR 

SUPPLY 


DISCHARGE  TO  ATMOSPHERE 


BOX  ENCLOSURE  AROUND  WICKET  GATE  ASSEMBLY  (ref  Fig.7) 


3.  STILLING  CHAMBER  DOWNSTREAM  OF  DRAFT  TUBE,  UPSTREAM  AIR  SUPPLY 


AIR 

SUPPLY 


VJ7 


DISCHARGE  TO 
ATMOSPHERE 


4.  STILLING  CHAMBER  AND  AIR  SUPPLY  UPSTREAM  OF  DRAFT  TUBE 


AIR 

SUPPLY 


DISCHARGE  TO 
ATMOSPHERE 


5,  STILLING  CHAMBER  WITH  SCREEN  DOWNSTREAM  OF  DRAFT  TUBE,  UPSTREAM  AIR  SUPPLY 


SCREEN 


AIR 

SUPPLY 


DISCHARGE 
TO  ATMOSPHERE 


- 5 ...  • . 
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TABLE  B2 


PRESSURE  PARAMETER, 


Figure  31.  The  effect  of  test  apparatus  geometry  on 
pressure  parameter  using  the  cylindrical 
draft  tube. 
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PRESSURE  PARAMETER 


Figure  32.  The  effect  of  test  apparatus  geometry  on 
pressure  parameter  using  the  elbow  type 
draft  tube. 
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Out  of  this  brief  study  comes  speculation  that  stilling  chambers 
placed  downstream  of  a turbine  runner  may  be  an  effective  means  of 
reducing  draft  tube  surge  in  hydroelectric  pump-turbines . Further  study 
is  required  in  this  area. 

A second  implication  has  been  suggested  by  H.  Falvey  of  The  Bureau 
of  Reclamation.  This  involves  the  geometry  of  model  test  loops  which 
have  stilling  chambers  located  up  or  downstream  of  model  pump- turbines . 
The  water  level  and  geometry  of  these  chambers  could  provide  surge 
characteristics,  based  on  model  tests,  that  are  inconsistent  with  those 
of  the  prototype. 
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